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Abstract. - We propose to realize the anisotropic triangular-lattice Bose- Hubbard model with 
positive tunneling matrix elements by using ultracold atoms in an optical lattice dressed by a fast 
lattice oscillation. This model exhibits frustrated antiferromagnetism at experimentally feasible 
temperatures; it interpolates between a classical rotor model for weak interaction, and a quantum 
spin-1/2 Xy-model in the limit of hard-core bosons. This allows to explore experimentally gapped 
spin liquid phases predicted recently [Schmied et al., New J. Phys. 10, 045017 (2008)]. 
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Introduction. — Frustrated quantum antiferromag- 
netism can give rise to extraordinarily rich physics pp. It 
is not only supposed to play a crucial role for the prop- 
erties of high-T c superconductors, but it is also an inter- 
esting subject on its own, and has potential applications 
in topological quantum information processing and stor- 
age [2J. Apart from the possibility of keeping a classical 
Nccl-ordered (staggered or spiral) spin configuration with 
long-range order, the spins of a quantum antifcrromagnet 
can also form singlet pairs ("valence bonds"), such that 
the spin-rotation symmetry is not broken. These singlets 
either order spatially to form a valence bond solid or the 
system's state is a superposition of many singlet coverings, 
neither breaking translations! nor spin-rotational symme- 
try. The latter is termed a resonating valence bond spin 
liquid (SL). A gapped SL with an exponential decay of spin 
correlations is expected to exhibit non-local topological or- 
der being immune against local perturbations. It also can 
possess anyonic excitations. This makes such topological 
SL states candidates for robust quantum memories and 
processors [2J. Alternatively, a critical SL is characterized 
by a huge density of low-lying excitations and a power-law 
decay of spin correlations. 

Since frustrated quantum antiferromagnets are hard to 
simulate (path-integral Monte-Carlo methods fail) and 
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clean solid-state realizations are not available, it is de- 
sirable to study these exotic many-body systems with ul- 
tracold atoms in optical lattice potentials [5] providing 
both clean conditions and far-reaching control. The most 
straightforward cold atom implementation of a quantum 
magnet is to create a Mott-insulator of fcrmions in two 
different internal states forming a pseudo spin. However, 
the necessary low temperatures (smaller than the weak 
superexchange spin coupling) have not yet been achieved. 
Also spinlcss fcrmions at filling 2/3 in a not yet realized 
trimerized Kagome lattice resemble a quantum magnet [1] . 

In this letter we propose a different strategy for the 
realization of a frustrated quantum system with ultra- 
cold atoms that — in contrast to the aforementioned ap- 
proaches — can be pursued in existing experimental se- 
tups, at temperatures already reached. Our idea is to 
consider spinless ultracold bosonic atoms in a triangular 
optical lattice and to induce frustration via a sign change 
of the matrix elements describing tunneling between ad- 
jacent potential minima. As will be shown, such a sign 
change can be achieved effectively by dressing the sys- 
tem with a high-frequency elliptical lattice acceleration. 
In the hard-core boson limit of strong repulsive interac- 
tion, the physics is then described by the antiferromag- 
netic spin-1/2 ATF-model on the triangular lattice. For 
certain regimes of anisotropic coupling this model is ex- 
pected to show gapped SL phases [5|- 
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This letter is organized as follows. We start with a dis- 
cussion of the frustrated positive-hopping Bosc-Hubbard 
model that describes the system to be realized experi- 
mentally. In order to sketch the expected phase diagram 
we combine (i) recently published numerical data (based 
on PEPS as well as exact simulations) [5] valid in the 
limit of strong interaction with (ii) results obtained by 
starting from the limit of weak interaction and system- 
atically including quantum fluctuations (beyond Bogoli- 
ubov). Then we discuss the experimental realization of 
the model, putting emphasis on how to change the sign of 
the tunneling matrix elements via a fast elliptical lattice 
acceleration. It follows a part devoted to the preparation 
of the frustrated model's ground state in the presence of 
a trapping potential. Finally, before giving a brief con- 
clusion, we discuss possible experimental signatures of the 
expected phases. 

Positive-hopping Bose-Hubbard model on a tri- 
angular lattice. — Consider a sample of ultracold 
bosonic atoms in a deep triangular optical lattice that is 
forced inertially by moving the lattice rapidly along an 
elliptical orbit. According to the following section, the 
time evolution of such a system has a simple description. 
Integrating out the fast oscillatory motion on the short 
time scale T = 2n/u) of the elliptical forcing, one finds the 
system's evolution on longer time scales governed by the 
time-independent effective Bose-Hubbard Hamiltonian 
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Here bi and fii are the bosonic annihilation and number 
operators for Wannier states localized at the minima of 
the triangular lattice potential. The first term comprises 
tunneling between adjacent sites i and j with — this is the 
crucial point — matrix elements J?? that are smoothly 
tunable from negative to positive values by variation of 
the forcing strength^ The on-site terms are character- 
ized by the positive interaction parameter U and the local 
chemical potential /ij = fx — Vi including the trapping po- 
tential Vi. We consider the anisotropic lattice shown in 
Fig. QJa) with the Jff equal to either J or J' = aJ (as- 
suming a > 0). 

The homogeneous model (fii = /i) interpolates between 
a classical rotor and a quantum spin model: For weak 
interaction U <C n\J\, with a mean filling of n particles 
per site, the supcrfluid (SF) ground state can (locally) be 
approximated by Yii exp^i&l) | vacuum) with discrete or- 
der parameter ipi = y / n7exp(i(/?i). A homogeneous density 
n, = n is favored and the local phases ifi play the role of 
classical rotors assuming a configuration ip^ = q ■ Ti de- 
scribed by the ordering vector q. Antifcrromagnetic cou- 
pling J > implies Neel ordered phases tpi as depicted 
in Fig. QJb-e). We call such a state a Neel SF. When 
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Fig. 1: (color online) (a) Anisotropic triangular lattice con- 
sidered, with primitive vectors ai = de x , = d[(l/2)e x + 
(v / 3/2)e H ], as well as a.3 = -01 + a-z- The tunneling matrix 
elements Jff take values J and J' = aJ (with a > 0) along 
the solid and dashed bonds, respectively, (b) Reciprocal lat- 
tice with b = (47r/v / 3)d _1 . The first Brillouin-zone, centered 
at p = 0, is shaded. Considering antiferromagnetic coupling 
J > 0, we have marked the ordering vector q describing a Neel 
SF in the limit of weak interaction: For a > Qo q lies on one of 
the x-shaped crosses (being equivalent modulo reciprocal lat- 
tice vectors). This corresponds to a staggered configuration of 
the local phase angles ifii on the rhombic lattice of J'-bonds 
[shown in (c) with the ipi visualized by pointers]. Lowering a, 
at a = ao the position of q splits in a continuous way into two 
non-equivalent possible positions that separate symmetrically 
along the arrows drawn in (b). The phases (ft assume a spiral 
pattern with two possible chiralities; subfigure (d) corresponds 
to the isotropic lattice (a = 1 < ao) with q lying on one of 
the corners of the first Brillouin zone. Finally, in the fD limit 
(a = 0) only q x has a defined value that is marked by the 
dashed lines in (b). The phase pattern is staggered along the 
ID chains of J-bonds [as sketched in (e)]. 



a exceeds a value ao, spiral continuously transforms into 
staggered Neel order. While ao equals 2 for U/J = 0, it 
slightly decreases with increasing interaction, cf. Fig.^a). 

In the opposite limit of strong interaction U 3> n\J\, 
there are only two energetically favored site occupations, 
rii = [n] = g (the largest integer smaller than n) and 
rii = g + 1. Associating them with "spin up" and "spin 
down" , respectively, gives the Bloch-sphere representation 
\*&i,ipi) = cos(^j/2)|g) i + sin(^/2)exp(i( ) s i )|3 + I) i at each 
site. Replacing (g + 1 )~ 1 b i by the spin lowering operator 
(erf — i<rf)/2, one arrives at the XF-model 
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1 In our convention (ij) denotes an oriented pair of neighboring 
sites, (ij) jt (ji). 



with hi = ^(JH~ Ug), J* Y = Siij??, and af, a\ , Sf 
being spin- 1/2 Pauli operators at site i. The ground state 
of Hxy cannot be a product state like nil^»>V«) wrtn 
definite local phases (fi anymore, since \di,tpi) cannot be 
an eigenstate of both af and a\ . Viewed from the Bose- 
Hubbard perspective, increasing intcrparticle repulsion in- 
creases the fluctuations of the local phases (ft. While for 
ferromagnetic coupling J < the classical phase config- 
uration is supposed to survive the presence of quantum 
fluctuations in the spin-I/2 limit U 3> n\ J\, for antiferro- 
magnetic coupling J > recent simulations suggest that 
(for Yli(<*f) = 0) classical Neel order is not necessarily 
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Fig. 2: (color online) Sketch of the phase diagram of the 
anisotropic positive-hopping Bose-Hubbard model on the tri- 
angular lattice for half-odd-integer filling. The parameter plane 
is spanned by interaction strength U/(nJ) and anisotropy ra- 
tio a = J'/J. The data (in green) for the spin-1/2 limit 
\U/(nJ) 3> 1] are taken from Ref. [5]- We assume the SL 
phases to survive at small finite values of J/U, since they 
are protected by a gap. The behaviour at small U / (nj) cor- 
responds to results obtained within a generalized Bogoliubov 
theory, cf. Fig. [5Ja). 



preserved [5]: Along the cv-axis different Neel phases are 
separated by gapped SL phases with exponentially decay- 
ing spin correlations. The results of Ref. [5] are displayed 
along the upper edge of the phase diagram shown in Fig. [21 
In order to gain further insight into the physics of the 
frustrated positive-hopping Bosc- Hubbard model (HJ, we 
start from the classical limit of weak interaction U <C n\J\ 
(assuming a homogeneous system) and include quantum 
fluctuations by using the generalized Bogoliubov approach 
introduced in Ref. H] . For filling already moderately larger 
than 1, we can replace b i ~ exp[i((^i + 5(f>i)]y/ni + Sfii, 
where Shi = 8h\ and Sifi ~ Sip\ describe quantum fluctu- 
ations to the local particle numbers rii and phases tpi, re- 
spectively, with [Shi,5ipi] ~ iSij. While {{5hi) 2 )i/n 2 <C 1 
can be assumed, the phase fluctuations {(5(f>i) 2 ) diverge in 
the ID limit (a = 0) (as well as at finite temperatures) 
where only quasi-long-range order is possible. However, 
one can still expect the fluctuation of the relative phases 
((dipt — Sipj) 2 ) between neighboring sites i and j to be 
small. Expanding the Hamiltonian (TIJ up to second order 
in Shi/rii and (Sifi — Sipj), it will be quadratic in terms of 
new bosonic operators = y/nl[8hi / (2n,i) + iSipi] and d], 
and can be diagonalized by a Bogoliubov transform (keep- 
ing (6hi) = 0). When computing, e.g., correlations (p\b-) 
between distant sites i and j, one cannot treat (dtfii — Sipj) 
as a small quantity, but has to use Wick's theorem to eval- 
uate expectation values of all powers of 8(pi [B]. We aug- 
ment this analysis by taking into account also the (Wick- 
decomposed) quartic corrections to the Hamiltonian when 
minimizing the ground-state energy with respect to both 
the Bogoliubov coefficients and the ordering vector 
This self-consistent above-Bogoliubov correction, that we 
include using a numeric iteration scheme, is necessary in 



order to explain a shift of aa with increasing interaction^ 
Assuming homogeneous filling rii = n as well as q y = 0, 
the method sketched in the above paragraph leads to the 
following results: With increasing interaction/quantum 
fluctuations, ao decreases, i.e. the cv-domain of rhombic- 
staggered Neel order grows [thick line in Fig. |3ja)] . This 
"order by disorder" phenomenon pQ is in accordance with 
the spin-1/2 results of Ref. [3] (cf. upper edge of Fig. [5]). 
In contrast, a finite ev-interval of staggered ID quasi-long- 
range Neel order, also predicted for the spin model, is not 
found. We have used these findings to draw the lower 
part of the phase diagram of Fig. [2] The quasiparticle 
dispersion relation is gaplcss and phonon-likc for quasimo- 
mentum wave numbers p around p = q. However, when- 
ever spiral order is found, it is symmetric with respect to 
(p — q) — > — (p — q) only in the limit of small \p — q\. In 
contrast, the zero temperature quasimomentum distribu- 
tion, having sharp peaks at p = q + reciprocal lattice vec- 
tors, possesses reflection symmetry with respect to p = q. 
Further results, for n = 3.5, are shown in Fig. [3fa). An 
estimate for the range of validity of the approximation is 
given by the dotted and the dashed line. Above them the 
relative phase fluctuations between neighboring sites sep- 
arated by di and ct2, respectively, exceed a value taken to 
be 7r/4. The fact that the dotted line does not approach 
zero in the limit of decoupled ID chains (a — > 0) indi- 
cates that in this limit the approximation still captures 
the physics in ai-dircction (along the chains). Moreover, 
the dip around a = «o can be interpreted as a precursor 
of the SL phase predicted in Ref. [5] (cf. Fig. [5]). Finally, 
the thin solid line marks the interaction where the con- 
densate fraction f c = lim^.^^^ \(bjbj)\/n is reduced to 
0.75. Again, a sharp dip at a = ao is a hint at a quantum 
disordered phase in the limit of large interaction. 

Proposal for an experimental realization. Hav- 
ing discussed the physics of the triangular-lattice positive- 
hopping Bose-Hubbard Hamiltonian, let us turn to the 
realization of the model with ultracold atoms in a deep 
optical lattice. The sign change of the tunneling matrix 
element, from negative to positive values, shall be induced 
by dressing the system with a fast time-periodic lattice ac- 
celeration. For hypcrcubic lattices, such a dynamical mod- 
ification of tunneling has been predicted theoretically not 
only for single [8], but also for many interacting particles 
[5]. Moreover, it has been observed experimentally with 
ultracold atoms both in the weakly interacting regime (via 
the expansion of a Bosc-Einstein condensate QH]), as well 
as in the strong-coupling regime, where it has been used 
to induce the quantum phase transition from a SF to a 
Mott insulator (MI) and back 9,11]. However, the linear 
driving scheme used in the work just mentioned, with the 



2 In the case of spiral order (ir < d\q x \ < 2ir with q y = 0), one 
finds two solutions, q and q' = —q, and has to choose one of them. 



3 Taking into account self-consistently the quartic terms does not 
lead to a spurious gap in the quasiparticle spectrum, as it is the case 
within the standard Bogoliubov treatment [7]. This gaplessness is, 
thus, a feature of the generalized Bogoliubov expansion [6] in terms 
of fluctuations in density and relative-phase. 
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Fig. 3: (color online) (a) Generalized Bogoliubov theory for 
n = 3.5. Values of U/(n\J\) at which: spiral changes to 
rhombic- staggered Neel order (thick line) , relative phase fluctu- 
ations reach tt/4 for sites separated by ai (dotted line) and ai 
(dashed line), the condensate fraction has dropped to 0.75 (thin 
solid line), (b) Boundaries of the MI phases with integer fill- 
ing n in the jj,/U-J/U-plane, both in 2nd-order strong-coupling 
(solid lines) and meanfield (dashed lines) approximation. As a 
consequence of frustration, the MI double-lobes are larger on 
the antiferromagnetic side ( J > 0) of the phase diagram. The 
grey bubbles between the MI regions, indicating the expected 
gapped SL phases at half-odd-integer filling, are just sketched. 



system being forced sinusoidally along a single direction 
(chosen to be diagonal with respect to all symmetry axes 
in the case of a square or a cubic lattice), is not suitable 
for the triangular lattice geometry. In order to be able to 
manipulate the system in a symmetric way with respect to 
the three non-orthogonal lattice directions a\, a,2, and 03 
[Fig. EJa)], here we propose to use elliptical forcing. This 
includes isotropic circular as well as linear forcing. 

The driving scheme to be considered can be realized 
inertially by moving the lattice along an elliptical orbit 
x(t) = Ax c cos(ijt)e c + Ax s sin(u;t)e s in space, with an- 
gular frequency oj, orthogonal unit vectors e c and e s , as 
well as amplitudes Ax c and Ax s . The resulting inertial 
force in the lattice frame of reference reads F(t) = —mx = 
F c cos(uit)e c + F s sin(a;£)e s where m is the boson mass and 
Fc/s = muj 2 Ax c / s . Choosing ui and F c / S small enough to 
exclude transitions from the lowest to higher Bloch bands, 
one can describe the system in the lattice frame of refer- 
ence by the driven Bosc-Hubbard model 



H{t) = J2 J iMh + ? E - + 5>t(*) - ^ 



(3) 

Here Jy < are the bare tunneling matrix elements and 
Vi{t) = —Ti ■ Fit) oscillating on-site energies. 

We assume that hio is large compared to the energy 
scales given by interaction (U), tunneling (n|Jy|), and 
trapping — with neighbors i and j), i.e. that the 
forcing is fast with respect to the time scales governing 
the undriven model. Under these conditions the time evo- 
lution of the driven system's state \ip(t)) will be to good 
approximation of the form 



The unitary operator 
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where 



Wi(t) = J dr«i(r)-| 



di' f dr«i(r), (6) 
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just describes a periodically time-dependent shift by —mx 
of the whole system in quasimomcntum, Wi = r, • mx. On 
top of this simple oscillatory motion on the short time scale 
T = 2tt/ui, the time evolution on longer times is governed 
by the effective time-independent Hamiltonian H c ff shown 
in Eq. ([IJ, namely 



\4> e fi{t)} = exp 



-H eS t)\M0)}. 



The dressed tunneling matrix elements are given by 

' Ki, 
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Here Jo is the zero-order Bessel function and Kij = 
(F c e c ■ Tij) 2 + (F s e s ■ t^-) 2 the amplitude of the poten- 
tial modulation between site i and j, where = — r.j. 
Thus, in the lattice frame, apart from the superimposed 
fast oscillation in quasimomentum, the system behaves as 
the one described by H e ff. When measuring the momen- 
tum distribution of the system in the laboratory frame 
by taking time-of-flight absorption images, one will en- 
counter the periodic quasimomcntum distribution of \ip e fi) 
at rest, being enveloped by the momentum distribution of 
the Wannier wave function oscillating like mx. 

The result presented in the preceding paragraph relies 
on the separation of time scales as well as on time av- 
eraging. We have obtained it within the framework of 
quantum Floquct theory |12] by generalizing the approach 
introduced in Rcfs. [MT5] to elliptical forcing. The deriva- 
tion is based on stationary degenerate-state perturbation 
theory on the level of an extended Hilbert space including 
time as a coordinate. Here we just give a simple argu- 
ment making the fleff-description plausible: Transform- 
ing = U'\il>) leads to the new Hamiltonian H' = 
WHU — ihU^dtU). Accordingly, H' is obtained from 
H by subracting the oscillating potential terms oc Vi(t) 



and replacing J^. 



J^exp(i[Wi - W 3 ]/h). Now the 



w U(t)\i; eS (t)). 



(4) 



rapidly oscillating phase factors in the tunneling terms of 
H' can approximately be taken into account on time av- 
erage, H'(t) -> i J^dtH'(t) = ffeff, giving |V) « l^eff). 

A 2D triangular optical lattice can be realized by super- 
imposing three laser beams, all polarized in z-direction, at 
an angle of 2-7r/3 in the xy-plane, while a standing light 
wave in z-direction is used to create a stack of effectively 
two-dimensional systems. A further beam in z-direction 
allows to modify the trapping potential in the xy-plane. 
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The lattice motion can be realized by varying the rela- 
tive frequencies of the beams by means of acousto-optical 
modulators. For the purposes described above, an orbit 
x(t) = Ax c cos(ujt)e c + Ax s sin(cdt)e s is required, with 
Ax s / C on the order of a lattice constant and u)/(2n) being 
a few kHz. Starting from an isotropic undriven lattice with 
bare tunneling matrix elements Jy = J < and choosing 
e c / s — e x / y , one obtains effective tunneling matrix ele- 
ments ([8]) distributed as depicted in Fig. [Ha). Namely Kij 
reads K = d\F c \ and K' = d^/F 2 + 3F 2 /2 along the solid 
and dashed bonds, respectively, giving J = JJq \KJ{t\iS)\ 
and J' = JJq(K' /(Huj)) according to Eq. ©. This allows 
for any value of the anisotropy parameter a = J' /J. 

We have already implemented a triangular optical lat- 
tice in the laboratory, loaded it with ultracold 87 Rb atoms, 
and observed the transition from a SF to a ML Also a con- 
trolled motion of the lattice has been achieved. 

State preparation and role of trapping poten- 
tial. — For elliptical forcing there are no instants in time 
where U(t) is equal to the identity (i.e. with x = 0) like it 
is the case for linear forcing (F s = 0) at integer tw/(2n). 
Thus, it is not possible to "map" the state of an ini- 
tially unforced system on \ipcfi) by suddenly switching on 
the forcing. However, one can smoothly switch on the 
drive. According to the adiabatic principle for quantum 
Floquet states [Tl], \ipefi) can follow adiabatically when 
H c g is modified by the forcing, starting from \tpes) = \t/>) 
in the undriven limit 0[T2]. Before passing from the 
ground state of the undriven system {Jf^ = Jij < 0) to 
the positive-hopping regime ( Jff > 0) in the presence of a 
trapping potential, the lattice should be tuned very deep, 
such that U ^> no| J\ with filling no in the trap center. The 
system will form MI regions |15j with an integer number 
g of particles (depending on Hi/U) localized at each site. 
Different MI regions will be separated only by tiny inter- 
mediate domains of non-integer filling. In the MI phases 
the state can follow smoothly through the sign-change of 
J when the lattice acceleration is ramped up continuously 
in a next step. Moreover, in a deep lattice unwanted in- 
terband transitions are strongly suppressed. When the 
desired strength of the forcing is reached, in the center of 
the trap the MI has to be melted. This can be achieved 
both by decreasing the lattice depth (without leaving the 
regime of strong correlation U ~ no | J\) and by tuning the 
chemical potential in the trap center. The latter can be 
achieved by varying the trap, such that atoms are pushed 
into or pulled out of the center. 

We have studied the MI-to-SF (J < 0) and Ml-to- 
Nccl SF (J > 0) transition in the triangular lattice the- 
oretically. In the parameter plane spanned by fi/U and 
J/U, a strong-coupling expansion as described in Rcf. |16| 
gives the upper and lower boundary, fX p /U and i^h/U, 
of the MI phase with integer filling n = g. One finds 
^ P /U = g - {g + l)r) - gcpV 2 + 0(t] 3 ) and f i h /U = (g - 
1 ) + QV + (.9 + 1 ) c hif + C ( ? 7 3 ) ■ The expansion parameter is 
given by r\ = —e{q)/U = w\J\/U with e{q) = —\J\w being 



the single-particle dispersion relation e(p) = 2 J[cos(dp x ) + 
2a cos(dp x /2) cos(y/3dp y /2)] evaluated at its minimum q. 
The coefficients read c p = g + 1 - (5g + 4)(1 + 2a 2 )/w 2 
and Ch = g— (5y-l- 1)(1 + 2a 2 )/w 2 . Here w directly reflects 
frustration; while w — 4a + 2 for ferromagnetic J < 0, it 
is smaller for antifcrromagnctic J > 0, namely w — a 2 + 2 
for < a < 2 and w = 4a — 2 for a > 2. As a conse- 
quence, the MI regions extend to larger values of \ J\/U 
on the frustrated side of the J — /i plane. This can also 
be observed in Fig. EJb) displaying the phase diagram for 
a = 1.3. Moreover, the transition from ID like concave 
phase boundaries c p /h < to square lattice like convex 
ones Cp/h > happens at noticeably larger a in the case 
of fustration. Namely it occurs for a between 0.03 and 0.13 
(2.3 and 4.2) when J < (for J > 0). For convex bound- 
aries, also fi p/h /U~ i{2 3 -l-77±[l-2(2(7 + l)r / + 77 2 ] 1 / 2 } 
obtained within meanfield approximation (cf. Refs. [15] ) 
can be expected to provide a reasonable description. 

The phase diagram plotted in Fig. [3£b) shows: the 
smaller n\J\/U gets, the smaller get the intervals of /i/U 
with non- integer filling (i.e. outside the MI lobes). In or- 
der to reach the strong coupling limit U 3> n\ J\ at non- 
integer filling g < n < g + 1 [where the spin- 1/2 descrip- 
tion (J2J with non-trivial polarization applies], the varia- 
tion of (i.e. of Vi) must be smaller than At|f +1 ' ) — ^p 9 ' ~ 
2(g + l)w\ J\ over an appreciable number of sites. Such a 
situation, where also the gapped SL phases are supposed 
to appear, can be achieved in the center of a shallow trap. 
Note that the presence of a (shallow) controllable trapping 
potential is definitely desirable: tuning its depth allows to 
manipulate the chemical potential/filling in the trap cen- 
ter. With respect to the chemical potential, the gapped 
SL phases, expected for large interaction and a near 0.5 
or 1.3 [cf. Fig. [5], would appear as incompressible regions 
at half-odd-integer filling. In Fig. EJb) we have sketched 
these phases (shaded in grey); they show up as "bubbles" 
on the frustrated side between the MI regions. 

Experimental signatures of frustration. Ex- 
perimental signatures of the Neel SF are sharp quasimo- 
mentum peaks at p = q + reciprocal lattice vectors [cf. 
Fig. [Ub)] . In the case of spiral order, the ordering vector q 
can take two different values; when measuring (or already 
before) the system will spontaneously choose one of them. 
For a whole stack of 2D systems, the measurement will 
average over both quasimomentum distributions, unless 
there remains a finite coupling between the 2D-layers es- 
tablishing the same order everywhere. Also the predicted 
downshift of the anisotropy ratio ao (where spiral continu- 
ously transforms into rhombic-staggered Neel order) with 
increasing interaction/lattice depth [cf. Figs. [5] and [3]^ a)] 
can be investigated experimentally The growth of the 
staggered a-domain with increasing quantum fluctuations 
is an example for "order by disorder" [T]. Another mea- 
surable consequence of frustration is the extension of the 
MI phases to larger values of \ J\/U [cf. Fig. Efb)] . Due to 
the lack of long-range order, the MI does not show sharp 
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peaks in the single- but rather in the two-particle momen- 
tum distribution (noise correlations) |T7] . This is also true 
for the gapped SL phases, being the most striking implica- 
tion of frustration expected. Thus, in order to distinguish- 
ing the SL from the MI experimentally, one should search 
for structures in the momentum distribution beyond sharp 
peaks. The SL can feature a pattern in the momentum dis- 
tribution on the scale of a Brillouin zone (i.e. the inverse 
lattice spacing n/d), reflecting derealization of particles 
on pairs of neighboring sites forming "singlets" . Apart 
from that, single-site resolved measurements [T5] clearly 
distinguish between MI and SL by number fluctuations. 

Conclusion and Outlook. — We have proposed to 
realize the positive-hopping Bosc-Hubbard model with a 
system of ultracold spinlcss atoms in a deep triangular op- 
tical lattice dressed by a rapid elliptical acceleration. Our 
scheme allows to experimentally investigate the physics 
of a frustrated quantum system under the clean and con- 
trolled conditions provided by ultracold atoms. Since frus- 
tration is induced to motional bosonic degrees of freedom, 
it is experimentally possible to reach temperatures that 
are low compared to the energy scales governing the sys- 
tem. The model smoothly approaches a quantum spin-1/2 
XY-model in the deep-lattice limit of strong interaction. 
In order to draw the phase diagrams shown in Figs. [5] 
and [2Kb) we have combined results from different ap- 
proaches: (i) numerical simulations applying to the spin- 
1/2 limit of strong interaction at half-odd-integer filling 
(published recently in Ref. [5]), (ii) an above-Bogoliubov 
theory valid in the limit of weak interaction, and (iii) ana- 
lytical strong-coupling as well as meanfield results for the 
limit of strong interaction at integer filling. Expected are 
supcrfluid phases showing staggered or spiral Neel order, 
Mott insulator phases having integer filling, and gapped 
spin-liquid phases at half-odd-intcger filling. We have also 
described how the positive-hopping regime can be reached 
adiabatically, if initially the system is prepared in the 
usual negative-hopping ground state. Finally, experimen- 
tal signatures of the different phases have been discussed. 
In conclusion, using an existing setup the experiment pro- 
posed here can provide novel information about a frus- 
trated quantum system. 

We have restricted our analysis to the triangular lat- 
tice geometry that we have implemented already in the 
laboratory. However, the route described here, namely 
(i) realizing a positive-hopping Bosc-Hubbard model with 
ultracold atoms by dressing the system by a fast ellipti- 
cal lattice acceleration and (ii) approaching the physics 
of a spin-1/2 XY-mode\ in the limit of strong interaction, 
applies equally to other two-dimensional non-bipartite lat- 
tices such as the Kagomc lattice. This opens perspectives 
for interesting future research. 
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